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ABSTRACT 

The Herschel OB young stellar objects survey (HOBYS) has observed the Rosette molecular cloud, providing an unprecedented view of its star 
formation activity These new far-infrared data reveal a population of compact young stellar objects whose physical properties we aim to charac- 
terise. We compiled a sample of protostars and their spectral energy distributions that covers the near-infrared to submillimetre wavelength range. 
These were used to constrain key properties in the protostellar evolution, bolometric luminosity, and envelope mass and to build an evolutionary 
diagram. Several clusters are distinguished including the cloud centre, the embedded clusters in the vicinity of luminous infrared sources, and the 
interaction region. The analysed protostellar population in Rosette ranges from 0.1 to about 15 M© with luminosities between 1 and 150 L , which 
extends the evolutionary diagram from low-mass protostars into the high-mass regime. Some sources lack counterparts at near- to mid-infrared 
wavelengths, indicating extreme youth. The central cluster and the Phelps & Lada 7 cluster appear less evolved than the remainder of the analysed 
protostellar population. For the central cluster, we find indications that about 25% of the protostars classified as Class I from near- to mid-infrared 
data are actually candidate Class objects. As a showcase for protostellar evolution, we analysed four protostars of low- to intermediate-mass in a 
single dense core, and they represent different evolutionary stages from Class to Class I. Their mid- to far-infrared spectral slopes flatten towards 
the Class I stage, and the 160 to 70yt/m flux ratio is greatest for the presumed Class source. This shows that the Herschel observations characterise 
the earliest stages of protostellar evolution in detail. 
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The HOBY S Hers chel (JPilbratt et al.l l2010h imaging survey 
(Motte et alJ I2010I) is a key programme for studying the 
sites of OB star formation within a distance of 3kpc. The 
parallel mode scan map ob servations, ca rried out with the 
PACS (iPoglitsch et al.l [20lob and SPIRE dGriffin et al.l l2010t 
Swi nyard et al.l l2010l) instruments, provide maps of massive 
cloud complexes in the wavelength range of 70 to 500 /mi, and 
thus a census of the different stages of star formation from 
prestellar cores to evolved young stellar objects (YSOs). As an 
excellent example of a massive cloud complex associated with 
- and under the influence of - an OB star cluster (NGC 2244), 
the Rosette molecular cloud was observed for HOBYS durin g 
the Herschel Science Demonstration Phase (JMotte et al.ll2010l) . 
The influence of NGC 2244 on th e cloud complex is di scussed 
in the accompanying pap er by [Schneide r et al.l ([2010), while 



iDi Francesco et al.l (120101) assesses the clump population. For 



consistency with Schneider et al. ( 20101) . the distance of Rosette 



* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im- 
portant participation from NASA. 

** Figures 5-10 are only available in electronic form via 
http://www.edpsciences.org 



adopted here is 1.6 kpc. Many previous studies have targeted 
the Rosette complex. Ongoing star formation across the molec- 
ular cloud is traced by the presence of several luminous IRAS 
sources t hat are associated with massive clumps seen in CO 
emission (ICox et al.lll990UWiliiams et"al . 1995; S chneider et al.l 
Il998l) . Towards these dense re gions, the embedded c lusters PL1 
to PL7 have been identified by Phel ps & Ladal (1 19971) from near- 
infrared o bservations. The l is t of clusters in Rosette wa s ex- 
tended by Ei & Smithl (120051) . iRoman-Zuniga et al.l (120081) . and 
Poul ton et al.l (120081) using near-infrared and Spitzer observa- 
tions. 

In this paper, we make use of the unprecedented spatial res- 
olution and sensitivity of Herschel to establish a sample of com- 
pact far-infrared sources that represent protostellar objects and 
to determine their fundamental properties (luminosities and en- 
velope masses) from their spectral energy distribution (SED). In 
the early phases of the collapse of a protostellar core and the 
initial accretion of matter onto a central protostar, the SED of 
a Class source is dominated by thermal emission from cold 
dust in the envelope (M env > M+. l Andre et al.ll2000h . The SED 
of a more evolved Class I source is shifted to mid-infrared wave- 
lengths, indicating a comparatively less massive, hotter envelope 
(Menv < M+). To investigate the evolutionary stage, the proto- 
stellar envelope mass is usually compared to the bolometric lu- 
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minosity, which serves as a proxy for stellar mass. The Herschel 
observations for the first time cover the peak of the protostellar 
SED thus constraining the evolutionary stage of early-to-evolved 
YSOs. This overcomes the previous difficulties distinguishing 
Class and Class I sources using near- to mid-infrared data. 

2. Observations 

The Rosette molecular cloud was observed by Herschel on 
October 20, 2009 in the parallel scan map mode (scanning speed 
of 207sec) simultaneously with SPIRE at 250/350/500 /mi and 
PACS at 70/1 60 /zm. Two perpendicular scans (consisting of 
parallel scanlegs interspersed by turn-arounds) were taken to 
cover a SPIRE/PACS common area of 1° x 1°. The data are re- 
duced with scripts developed in HIPEQ (Ott 2010, version 2.0 for 
SPIRE and version 3.0 for PACS). The PACS data are deglitched 
from cosmic ray impacts with the HIPE second-level method and 
then high-pass filtered with a scanleg filter width to preserve the 
extended emission up to the map size scale. The combined scans 
are finally projected using the HIPE MadMAP implementation 
with the noise table Invntt version 1. For details of the SPIRE 
data reduction, see ISchneider et alJ (120101) . The maps are flux- 
calib rated accordi ng to the corr ection factors of ISwinyard et al] 
d2010b and Poglit sch et al] d2010h . compared to 2MASS to cor- 
rect for a ~6" pointing offset and a systematic offset of ~4" 
between SPIRE and PACS. The entire set of Herschel maps is 
shown in lMotte et al] (120101) . and Fig. |4l (available online) shows 
the 70 and 160 /mi maps. 




6 h 34 m 30 s 



20 s 
R.A. (J2000) 



Fig. 1. PACS 70yum of the Rosette molecular cloud centre. The 
region harbours the embedded clusters PL4 (northwest), PL5 
(southeast), and a concentration of compact Herschel sources. 



3. Results and analysis 

Rosette harbours several luminous far-infrared sources that rep- 
resent candid ate high-ma ss protostars with AFGL 961 being the 
brightest (see Mott e et alj|2010l) . In their vicinity, the Herschel 
70 and 160 /mi maps reveal a population of compact sources that 
likely represent YSOs of low to intermediate mass. The most 
prominent cluster is the one towards the Rosette molecular cloud 
centre (see Fig.Q]). The emission traces heated protostellar enve- 
lope material; in contrast, starless cores are not expected to be 
detected as compact sources at 70 /mi. That Herschel resolves 
individual protostars, even though the Rosette region lies at an 
intermediate distance, allows us to study the physical properties 
of the protostar population based on unprecedented far-infrared 
data. 



3.1. Identification of compact protostellar sources 

To identify the positions of the compact Herschel pro- 
tostars, we made use of the 70jim MRE-GCL catalogue 
( Motte et all I201 Oh derived using the method of Motte et all 
(2007). Sp atial scales larger than 0.5 pc were filtered out us- 
ing MRE (Starck & Murtagh 2006), and the Gaussclumps pro- 
gramme applied (Kram eret al] 1 1998k . We focused on the pro- 
tostellar objects clearly detected at 70 /mi with FWHM < 15" 
and derived the 70 and 160 /mi flux measurements using aperture 
photometry. Sources that are faint or that appear non-singular, 
either of which prevents a useful measurement, were excluded. 
The aperture diameters were chosen to correspond to a phys- 
ical scale of 0.1 pc (corresponding to a source FWHM size of 
roughly 0.05 pc) and the background is determined in adjacent 



1 HIPE is a joint development by the Herschel Science Ground 
Segment Consortium, consisting of ESA, the NASA Herschel Science 
Center, and the HIFI, PACS, and SPIRE consortia. 



annuli. Figures Q] and [3] show that protostars are resolved at 70 
and 160/zm but blend together in the submillimetre, preventing 
a dust temperature measurement on the protostar scale. Thus we 
adopted temperatures determined by Motte et al . (2010) derived 
from the integrated emission on larger scales, using greybody fits 
with a dust emissivity index of /? = 2 and assuming Tj = 20 K 
for sources where no estimate is available. This may introduce a 
bias towards low temperatures as the protostellar envelopes are 
expected to be warmer than their surroundings. The envelope 
masses were estimated by scaling the greybody curve (of dense 
cores) to the 160 /mi fluxes (of protostellar envelopes). Eighty 
eight protostars are the basis for the further study described be- 
low. Due to the compactness criterion, this sample does not in- 
clude the high-luminosity, high-mass protostellar objects (e.g. 
AFGL 961). 

Partly based o n Spitzer IRAC and MIPS data of the same 
region obtained bv lPoulton et al.l (120081) and on 2MASS, a cat- 
alogue supplied by Balog et al. (in prep., hereafter referred to 
as Spitzer catalogue) was used to search for source counterparts 
in the near- and mid-infrared. Considering -1900 sources cata- 
logued at 24//m, 17 Herschel sources have no counterpart. The 
inspection of the maps shows that most of them are either not 
covered by the MIPS observations or are not included in the cat- 
alogue owing to artifacts in the Spitzer maps. Interestingly, 5 
sources have only been detected with Herschel but not at 24 /mi 
(completeness limit -0.5 mJy). The Spitzer catalogue also in- 
cludes a classification of source type using the IRAC colours 
Gutermu th et al] (120081) . The compiled flux measurements that 
include the Spitzer catalogue are used to estimate bolometric lu- 
minosities by integration over the SEDs, where we use the sub- 
millimetre fluxes of a greybody with the adopted dust temper- 
ature. The current lack of resolved submillimetre flux measure- 
ments means that the mass derivation relies on assumptions of 
the dust temperature and emissivity, and we estimate that the 
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Table 1. Protostar subsamples seen by Herschel in Rosette 



O Central cluster 

O PL7 

O other Rosette protostars 



10" 



Fig. 2. Envelope mass versus bolometric luminosity diagram 
for the sample of Herschel protostars in Rosette. Evolutionary 
tracks for stellar masses between 0.2 and 20 M© are included. 
The solid line corresponds to 50% of the mass accreted; the 
dashed lines account for the estimated uncertainties of M env 
and Lboi- Open symbols: candidate Class protostars with 
La> 350//m/^bol > 1%; filled symbols: Class I protostars with 
Ai>350/im/A>oi ^ Ir- 



relative accuracy is roughly a factor of 2, and similarly for the 
luminosities. This will be improved by forthcoming studies. 

3.2. Evolutionary stage of the Herschel protostars in Rosette 

The derived envelope masses and bolometric luminosities are 
plotted in an evolutionary diagram shown in Fig. [2l For clar- 
ity, we omit the aforementioned error estimates. Such a dia- 
gram is proposed to trace the evolution of embedded protostars 
(Bontemps et al] 11996k ISaraceno et al.l 1199 6) because Class 
sources display relatively high envelope masses that decrease 
during the mass accretion, while the luminosities (protostellar 
core and accretion) increase significantly. E volutionary tracks 
are d isplayed fo r stars of different masses (cf. lAndre et al.l l2008: 
Bontem ps et al.ll2010h . The Rosette protostar sample occupies 
the low- to high-mass regimes in the diagram. C ompared to find- 
ings i n the nearby Aquila star-forming region dBontemps et al.l 
l2010h . the Rosette protostars are comprised of higher masses. 
The diagonal lines in Fig. [2] (see caption) indicate an approxi- 
mate border zone between envelope-dominated Class and star- 
domina ted Class I objects based o n the comparison of M Qm to 
M* (cf. lAndre & MontmerldFl994l) . 

A surprisingly large fraction of our sample (-2/3) falls into 
the candidate Class regime above these lines. A practical cri- 
terion i nferred from this d iagram is L^ > 35o^m/^bol > 1% f° r 
Class (I Andre et al .1120001) . In Fig. [2] we apply this criterion ob- 
tained from the previously compiled SEDs. It results in more 
intermediate-mass objects being classified as Class I, i.e., a more 
conservative Class assignment that we hereafter refer to as 
"candidate Class 0". The classification of objects lying near the 
border zone remains tentative because of the uncertainty in con- 
verting the M e nv/^* ratio into measurable quantities. In particu- 
lar, the Herschel photometry of nearby protostars will reduce this 
uncertainty, after data from several regions have been analysed. 
Until a detailed analysis of the completeness is available, our 



Region 



Associated Clusters # protostars 



Centre 


PL4/5, REFL08, E 


27 


Shell 


PL1,A 


5 


Monoceros Ridge/ 






Extended Ridge 


PL2,C 


11 


PL3 


- 


4 


AFGL961 


PL6 


6 


PL7 


G 


11 


Distributed 


- 


24 



Notes. Cluster names are PL for Phelps & Lada (1997), REFL for 
Roman-Zuniga et al. (2008), A. . . G for Poulton et al. (2008). 



finding thus remains preliminary. Nevertheless, it indicates that 
Herschel allows us to significantly extend the sample of known 
Class objects. 

Seven protostar subsamples are es tablished according t o 
their location (field denominations in ISchneider et al] 120101) . 
which are listed in Table [T] The distributed protostars con- 
tains sources located towards the tip of pillar structures 
(ISchneider et al] 120101) . Each subsample spans a wide mass 
range from 0.1 to about 10M o . In Fig. [2] the central cluster and 
PL7 subsamples are emphasised by coloured symbols. In the 2 
to 10M© range, the central cluster harbours a significant num- 
ber of sources in the Class regime with 4 to 30 L©. Also the 
PL7 cluster contains 3 sources with relatively low luminosities. 
This indicates that both the central cluster and the PL7 cluster 
are younger compared to the remaining protostars. 

3.3. The classification of protostars in the central cluster 

Among the 27 Herschel protostars in our central cluster sam- 
ple, the Spitzer classification (see Sect. 13.11) gives 12 Class I 
objects and one for Class II. For PL7, 5 out of 11 sources 
are classified as Class I. These classifications are added to the 
symbols in Fig. [2] as T and 'IF, for comparison to their lo- 
cation in the evolutionary diagram and the distinction by the 
La > 350/zm/^bol ratio. Notably, about half of the sources seen as 
Class I by Spitzer in these two subsamples correspond to sources 
with Lx > 350//m/^bol > 1% and thus represent candidate Class 
protostars. They are intermixed with the unclassified sources in 
the diagram that lack detections in one or several bands. This 
suggests that the classification based on the near- to mid-infrared 
data alone has to be partly revised in light of the Herschel mea- 
surement of the protostellar envelope. 

To overcome incompleteness in our Herschel catalogue and 
to provide a first census of how many sources classified using 
near- to mid-infrared data are affected, we have focused on the 
central cluster. We redefined the cluster area as a rectangular 
field towards the Rosette molecular cloud centre and selected 
a total of 86 catalogued 24 /mi sources in that field. Three 24 /mi 
sources were classified as probably extragalactic and are ex- 
cluded. We then inspected the Herschel 70 /mi map for coun- 
terparts. The field and the resulting detection statistics are given 
in Table [2] The chosen field excludes PL4 and its vicinity be- 
cause the extended emission there makes source identification 
very uncertain. Still there is a bias towards bright objects, and we 
estimate an uncertainty of three sources in total and one source 
per subsample. Based on visual inspection, we find a 70 /mi de- 
tection rate of about 50% for YSOs seen at 24 /mi. Due to their 
rising SED, we expect more 70 /mi detections for Class I than for 
Class II, and find about 3/4 compared to about 1/4. In the second 
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Table 2. Classification of protostars in the Rosette central cluster 





# total 


# Class II 


# Class I 


# unclass. 


24 //m sources 
Visible in 70yt/m 


83 

40 (±3) 


39 
10 (±1) 


26 

19 (±1) 


18 
11 (±1) 


In Herschel sample 
Candidate Class 


22 
14 


1 



12 

7 


9 

7 



Notes. The cluster area is defined here by 98.5° < RA < 98.6° 
and 4.25° < DEC < 4.4°. The table head classifications base on 
near- to mid-infrared measurements. Class candidates exhibit 

L ^>350//m/ L bol > 1% - 



b&2Qm Htmn 




1Q 13 



- Source 1 (Class I) 
Source 2 

- Source 3 

- Source 4 (24yu.m upper limit) 




Wavelength (/xm) 

Fig. 3. Four protostars in the central cluster at 24, 70, 160, and 
250 yum and their spectral energy distributions. 



row of Table [21 we list the corresponding numbers for the anal- 
ysed Herschel sample and the classification for the same field. 
Roughly, we include about half of the visible 70jim sources 
in the Herschel list. About 2/3 of these are Class candidates. 
This applies to 7 out of 18 unclassified sources, which indicates 
that many of the latter are Class protostars. About 25% of 26 
previously classified Class I sources are also Class candidates, 
possibly more. Notably, our statistical basis is low and will be 
improved in forthcoming studies. 

For illustration, we consider a particularly interesting large 
dense core revealed by Herschel which is resolved into several 
proto stars at 70/160/mi (Fig. [3]). For the whole core jMotte et al.l 
(2010) derive a dust temperature of 15 K. We focus on 4 sources 
clearly detected at 70/mi labelled 1 to 4 from west to east. The 
projected separation between two neighbours is around 0.14pc. 
At 24 /mi, only sources 1 to 3 are detected with decreasing 
brightness. Based on the Spitzer data, source 1 is classified 
as Class I, while the others remain unclassified owing to non- 
detection in one or more of the IRAC bands (sources 2 and 3) or 
non-detection in all near- to mid-infrared bands (source 4). At 70 



and 160yum source 4 shows up and the adjacent source 3 is the 
brightest of the group. The SED slopes (Fig. [3]) between 8 and 
70yum for sources 1, 2, and 3 become increasingly steep. Beyond 
70yum they are shallower, and the 160 to 70/mi flux ratio is the 
highest for source 4. The mid- to far-infrared colours allow us to 
distinguish the source evolutionary stages. 

The derived luminosities (M^ m pc , L^ > 350/mi/^bol) are 
15 L Q (2M , 1%) for source 1, 8L (3M Q , 3%) for source 2, 
8 L (5 M , 4%) for source 3, and 5 L Q (4 M Q , 6%) for source 4. 
We note that the assumption of a single dust temperature means 
the relative mass differences are not well constrained. Relatively 
similar in mass, the sources represent successive evolutionary 
stages from early Class (source 4) to "flat- spectrum" Class I 
(source 1). This is a first example of how Herschel resolves the 
early evolution of protostars in combination with Spitzer mea- 
surements. 



4. Implications 

This initial study of the protostellar population in Rosette us- 
ing the HOBYS observations already shows that Herschel pro- 
vides detailed insight into the earliest stages of the protostel- 
lar evolution at low to high masses, and in particular identifies 
thus far elusive Class objects. The now available Herschel pho- 
tometry will be the basis for establishing firm criteria (mid- to 
far-infrared colours) for distinguishing the evolutionary stages, 
including the observations of nearby regions where individual 
(low-mass) protostars are resolved in the submillimetre. Based 
on a complete protostar catalogue of the Rosette complex, we 
will examine the evolution of the star formation activity over the 
whole cloud also with respect to triggering (cf . ISchneider et al.l 
I201QI) . 
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Fig. 4. HOBYS Herschel 70 and 160 /mi maps of the Rosette molecular cloud (flux unit: MJy/sr). See lMotte et al.l (l201Qh for details 
on the sensitivity of the entire Herschel data set. 



